Discovery of X-Ray Burst Oscillations from a Neutron-Star 
X-Ray Transient in the Globular Cluster NGC 6440 
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ABSTRACT 

We report the discovery of millisecond oscillations in an X-ray burst from the X-ray transient 
SAX J1748. 9-2021 in the globular cluster NGC 6440. Oscillations at a frequency of 409.7 ±0.3 Hz 
were present in one of nine X-ray bursts observed with the Proportional Counter Array on the 
Rossi X-ray Timing Explorer during the outburst which occurred in 2001. The burst was one of 
the two dimmest and had the longest duration and decay time. The average peak luminosity of 
two bursts showing radius expansion is (3.6 ± 0.4) x 10 38 ergs _1 , consistent with the Eddington 
luminosity for a IAMq and 10 km radius neutron star burning hydrogen-poor matter. We 
speculate that the dichotomy observed between sources with burst oscillations at once versus 
twice the frequency difference of kHz quasiperiodic oscillations in the persistent emission may be 
related to the magnetic field geometry of the neutron stars. 



Subject headings: Galaxy: globular clusters: individual (NGC 6440) 
- stars: individual (SAX J1748.9-2021) 



stars: neutron — X-rays: binaries 



1. Introduction 

The discovery of millisecond oscillations in ther- 
monuclear X-ray bursts provided the first direct 
measurements of the spin frequencies of neu- 
tron stars in low-mass X-ray binaries (LMXBs) 
(Strohmayer et al. 1996; Strohmayer & Markwardt 
2002). The presence of neutron stars rotating with 
millisecond periods in NS-LMXBs was suggested 
more than a decade previously (Alpar et al. 1982) 
and is key in understanding the relation of LMXBs 
and millisecond radio pulsars. 

X-ray burst oscillation sources form two dis- 
tinct classes: "fast oscillators" with frequencies 



close to twice the frequency difference of the kHz 
quasiperiodic oscillations (QPOs) seen in the per- 
sistent emission, and "slow oscillators" with fre- 
quencies near the difference of the kHz QPO fre- 
quencies (White & Zhang 1997). If the kHz QPO 
frequency difference is close to the spin frequency 
of the neutron star, then the slow oscillators have 
a burst oscillation frequency close to the spin 
frequency, while the fast oscillators produce sig- 
nals at twice the spin frequency. The fast os- 
cillators produce burst oscillations predominately, 
but not exclusively, in photospheric radius expan- 
sion bursts, while the slow oscillators produce os- 
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dilations in bursts both with and without pho- 
tospheric radius expansion (Muno et al. 2001). 
This indicates that the difference between the two 
classes is a physical difference and not a selection 
effect. The nature of this physical difference is 
unknown. 

If the fast oscillators indicate twice the spin fre- 
quency and the slow oscillators indicate the spin 
frequency, then the spin frequency distribution of 
neutron stars in X-ray bursters is clustered with 
a maximum value near 400 Hz. Clustering in the 
spin frequency distribution has been variously in- 
terpreted as evidence for very efficient mass ejec- 
tion by rapidly rotating neutron stars in bina- 
ries (Burderi et al. 2001), a gravitational radia- 
tion limit on the spin-up of accreting neutron stars 
(Bildstcn 1998) and for the existence of a phase 
change in superdense nuclear matter (Glendenning 
& Weber 2001). 

Increasing the sample of sources from which 
burst oscillations are detected is crucial for deter- 
mining the physical nature of the difference be- 
tween the fast and slow oscillators and in accu- 
rately measuring the spin frequency distribution 
of neutron stars in LMXBs. Here, we describe ob- 
servations made with the Rossi X-Ray Timing Ex- 
plorer (RXTE; Bradt, Rothschild, & Swank 1993) 
following the detection of X-ray bursts from the 
transient source SAX J1748. 9-2021 located in the 
globular cluster NGC 6440. We report the discov- 
ery of millisecond oscillations in one X-ray burst. 
Because the source is located in a globular cluster, 
the distance is known and accurate luminosities 
can be calculated. We find that the peak lumi- 
nosity of the radius expansion bursts is consistent 
with the Eddington limit for a 1.4M Q and 10 km 
radius neutron star with hydrogen-poor matter. 
We describe the source and our observations in 
§2, the X-ray bursts in §3, the persistent emission 
in §4, and conclude in §5 with a few comments on 
the nature of the difference between the fast and 
slow X-ray burst oscillators. 

2. Observations of NGC 6440 

NGC 6440 is a globular cluster at a distance of 
8.5 ±0.4 kpc (Ortolani, Barbuy, & Bica 1994) and 
located near the Galactic center. Bright, transient 
X-ray emission coincident with NGC 6440 was first 
detected with OSO-7 and UHURU in December 
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Fig. 1 — RXTE/ ASM light curve of the 2001 out- 
burst from NGC 6440. The points are one day 
averages of the ASM count rates. The diamonds 
indicate the times of the PCA observations. 

1971 (Markert et al. 1975; Forman, Jones, Tanan- 
baum 1976). Strong X-ray outbursts were found, 
again, in August 1998 (in 't Zand et al. 1999) and 
August 2001 (in 't Zand ct al. 2001). The peak lu- 
minosities of the outbursts were near 10 37 ergs _1 . 
The durations of the outbursts were from one to 
a few months. The 1998 and 2001 outbursts are 
from the same object (in 't Zand ct al. 2001) iden- 
tified with the optical star referred to as V2 in 
(Vcrbunt et al. 2000). Whether the 1971 outburst 
is from the same object or one of the other, cur- 
rently quiescent, LMXBs in NGC 6440 (Pooley 
et al. 2002) is unknown. Type-I X-ray bursts de- 
tected in 1998 and 2001 with the Wide-Field Cam- 
era (WFC) on BeppoSAX and a high resolution 
image of the cluster obtained with Chandra while 
the source was active during 2001 clearly identify 
the source SAX J1748.9-2021 with V2 and as an 
accreting neutron star (in 't Zand et al. 2001). 

The RXTE All-Sky Monitor (ASM) light curve 
of the 2001 outburst is shown in Fig. 1. The out- 
burst started around MJD 52129 (8 August 2001), 
as reported by the MIT ASM team. This initial 
event was relatively weak at its peak and lasted 
only 12 days. After a pause of about 20 days, 
NGC 6440 again became bright in x-rays. The 
x-ray emission rose gradually over about 10 days 
and then jumped rapidly to the peak flux level 
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Fig. 2. — Dynamical power spectrum of burst 9 
with the burst light curve (count rate) superim- 
posed. The contours are at Leahy powers of 12, 20, 
and 28. The contours are generated from power 
spectra for overlapping 3 s intervals of data with 
the power plotted at the mid-point of the interval. 
The left vertical axis gives frequency values for the 
contours. The right axis gives count rates, for the 
sum of the 3 PCUs on during the burst, for the 
light curve. 



near 12 ASM c/s, which corresponds to a flux near 
3.2 x 10~ 9 crgcm~ 2 s _1 in the 2-10 keV band or 
a luminosity of 3 x 10 37 ergs~ 1 for a source lo- 
cated at the distance of NGC 6440. High flux 
levels persisted for about 20 days, and then the 
source rapidly declined and became undetectable 
with the ASM. 

The X-ray bursts found with the WFC during 
the 2001 outburst triggered an RXTE Targct-of- 
Opportunity program which led to an observation 
on 2001 September 28 and multiple observations 
during October 7-16. The total observation time 
was 81 ks. Our observations occurred during the 
bright phase of the outburst and in the initial part 
of the final decline. Data were obtained from the 
Proportional Counter Array (PCA) in the same 
spectral and timing modes described in Kaaret et 
al. (2002). 



3. X-Ray Bursts 

We used the Standard 1 data which has 0.125 s 
time resolution and no energy information to 
search for X-ray bursts and found nine bursts, 
see Table 1. We extracted energy spectra for the 
bursts from event mode data with 64 channels of 
energy resolution using all Proportional Counter 
Units (PCUs) on during each burst and all lay- 
ers. We include data from PCU0 because the 
problems with background rejection and modeling 
caused by the loss of its propane layer are not 
important for spectral analysis of these relatively 
bright bursts. We note that we did use the most 
recent response matrices, which include a correc- 
tion for the spectral response due to the loss of 
the propane layer. Spectra were accumulated for 
0.25 s intervals and corrected for deadtime effects. 
A spectrum from 10 s of data preceding each burst 
was subtracted to eliminate the contribution of the 
persistent emission. 

We fitted the resulting spectra in the 3-20 keV 
band with an absorbed blackbody model with the 
column density fixed to 5.9 x 10 21 cm~ 2 , the value 
equivalent to the optical reddening (Pooley et al. 
2002). The results are shown in Table 1. The 
bolometric flux was calculated from the spectral 
fit. The equivalent blackbody radius quoted in 
the table was calculated assuming a distance of 
8.5 kpc and the uncertainty docs not reflect the 
distance uncertainty. We note that these charac- 
terizations are only approximate since significant 
deviations from blackbody spectra have been de- 
tected from some sources and reprocessing of the 
radiation in the neutron star atmosphere alters the 
observed temperature and inferred radius from the 
true values (Lewin et al. 1995). 

Bursts 1 and 5 show evidence for radius ex- 
pansion including a sharp increase in radius of 
more than 20 km and a simultaneous decrease 
in temperature of more than 0.7 keV. The peak 
fluxes for these two bursts are consistent within 
the errors and their average peak luminosity is 
(4.3 ± 0.5) x 10 38 ergs _1 where the uncertainty 
includes the distance uncertainty. The fluxes re- 
ported from the RXTE PCA with PCUs 0, 1, 
2, and 3 operating (which was the case during 
bursts 1 and 5) are systematically a factor of 1.2 
higher than from other instruments. Assuming 
that this factor represents an error in the abso- 
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Time 


Peak Flux 
(10~ 8 crg cm~ 2 s _1 ) 


Peak Radius 
(km) 


Decay Time 

(8) 


Duration 

(s) 


1 


Oct 8 at 09:20:50 


5.1 ±0.5 


48.3 ± 1.2 


3.7 


11 


2 


Oct 8 at 11:15:05 


3.8 ±0.5 


14.3 ±0.9 


3.9 


12 


3 


Oct 10 at 06:25:36 


4.6 ±0.5 


19.9 ± 1.3 


4.3 


12 


4 


Oct 10 at 09:25:03 


2.5 ±0.4 


9.1 ±1.1 


9.4 


24 


5 


Oct 11 at 07:49:37 


4.8 ±0.5 


23.5 ±1.5 


3.9 


11 


6 


Oct 11 at 09:40:37 


3.6 ±0.5 


12.7 ±0.9 


6.6 


18 


7 


Oct 13 at 07:43:35 


3.7 ±0.5 


18.1 ± 1.4 


4.4 


14 


8 


Oct 16 at 04:58:12 


4.0 ±0.6 


9.3 ±0.7 


6.8 


21 


9 


Oct 16 at 08:40:28 


2.5 ±0.4 


8.5 ±0.9 


13.0 


36 



Table 1: Properties of X-Ray Bursts. The table gives for each burst: the time (UTC) in the year 2001 at the 
start the burst, the bolometric peak flux, the maximum blackbody radius, the 1/e-folding time in the tail of 
the decay, and the duration defined as the interval over which the flux is greater than 10% of the peak flux. 



lute calibration of the PCA (Tomsick et al. 1999), 
we reduce the flux by this factor and find an 
average peak luminosity for these two bursts of 
(3.6 ± 0.4) x 10 38 ergs _1 which is consistent with 
the Eddington luminosity for a 1.4M and 10 km 
radius neutron star with hydrogen-poor matter 
(Kuulkers et al. 2002). The other bursts have 
lower peak fluxes. Burst 3 has a maximum ra- 
dius close to 20 km and exhibits a temperature 
drop of 0.8 keV. Bursts 2, 6 and 7 show weak evi- 
dence of radius expansion with temperature drops 
of 0.4-0.7 keV but radius expansion of less than 
20 km. Bursts 4, 8 and 9 show no evidence of ra- 
dius expansion, no temperature drops, and decay 
more slowly than the other bursts. 

Using high time resolution data (merged event 
lists from the event and burst catcher modes) with 
no energy selection, we computed power spec- 
tra for overlapping 3 s intervals of data, with 
0.25 s between the starts of successive intervals, 
and searched for excess power in the range 100- 
1000 Hz. We found oscillations in burst 9, with 
a maximum Leahy normalized power of 38.7 at a 
frequency of 409.7 Hz, occurring in the burst de- 
cay about 5 s after the burst peak. Given 2700 
independent frequencies in the 900 Hz interval 
searched, the signal power corresponds to a proba- 
bility of chance occurrence of 1.1 x 10~ 5 , equivalent 
to a 4.4(7 detection. The dynamical power spec- 
trum is shown in Fig. 2. The peak rms amplitude 
of the oscillation is near 6%. Burst 9 is one of the 
two dimmest bursts and has the longest duration 
and exponential decay time. 



4. Persistent emission 

We analyzed the persistent emission following 
the procedures described in Kaaret et al. (2002). 
We used only data from PCU2 because this PCU 
was on during all of the observations. PCU0 was 
also on during all observations, but a fully reli- 
able background model for PCU0 subsequent to 
the loss of its propane layer is not available. The 
PCA light curve shows a high level of variability, 
up to 80% within one day. An X-ray color-color 
diagram is shown in Fig. 3 with each point repre- 
senting a 256 s integration. The points lie mainly 
in a single cluster, except for a few with low soft 
colors. A power spectrum for all of the data with 
soft color above 1.5 is shown in Fig. 4. The power 
spectrum has very low frequency noise and a so- 
called high frequency noise component as typically 
seen from atoll sources in the lower banana state 
(van der Klis 1995). The solid line superimposed 
on the data is a fitted model consisting of the sum 
of a low-frequency powerlaw, v 1 with 7 = —2.24, 
for the very low frequency noise and an exponen- 
tially cutoff powerlaw, ifte~ v / Vcnt w ^h ^ = q g2 
and v cut = 14.5 Hz. The RMS fraction is 5.7% 
in the 0.1-100 Hz range. The points with soft 
color below 1.5 have the lowest count rates and are 
from a single, contiguous time interval. A power 
spectrum of these data does not appear to have 
the very low frequency noise component, which 
would be consistent with identification as the is- 
land state. The RMS fraction is about 5% in the 
0.1-100 Hz range. 
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Fig. 3. — Color-color diagram of the persistent 
emission. Each point represents 256 s of data. The 
hard color is defined as the count rate in the 9.7- 
16 keV band divided that in the 6.4-9.7 keV band; 
the soft color is the 4.0-6.4 keV band rate divided 
by the 2.6-4.0 keV band rate. 

We searched for high frequency QPOs in each 
uninterrupted RXTE observation window and in 
combinations of the various data segments, includ- 
ing the sum of the island state data. We calculated 
averages of 2 s power spectra for all PCA events 
(2-60 keV) and for events in the 4.7-20.8 keV en- 
ergy band. We included events from all PCUs on 
during each observation. We did not find any sta- 
tistically significant signals. 

5. Discussion 

Our discovery of burst oscillations from the X- 
ray transient in NGC 6440 adds another source to 
the 10 sources previously known to exhibit burst 
oscillations. As described above, X-ray burst oscil- 
lation sources form two classes: "slow oscillators" 
versus "fast oscillators" with the distinction being 
whether the burst oscillation frequency is close to 
once or twice the frequency difference of the kHz 
QPOs in the persistent emission (White & Zhang 
1997). The phenomenology of the fast and slow 
oscillators is different with the fast oscillators be- 
ing much more likely to produce burst oscillations 




1.0 10.0 
Frequency (Hz) 



Fig. 4. — Low frequency power spectrum for all 
observations with soft color above 1.5. The power 
is RMS normalized. The solid line is a fitted curve 
described in the text. 

in photospheric radius expansion bursts which in- 
dicates that there is a physical difference between 
the two classes (Muno et al. 2001). If the kHz 
QPO frequency difference is close to the spin fre- 
quency of the neutron star, then the slow oscilla- 
tors produce signals at close to the spin frequency 
and the fast oscillators at twice the spin frequency. 

The fast oscillators have frequencies in the 
range 521-620 Hz, while the slow oscillators have 
frequencies in the range 270-363 Hz. The 410 Hz 
burst oscillation frequency of SAX J1748.9-2021 
lies between the two classes. The detection of 
oscillations during a non-radius expansion burst 
and the absence of detections in radius expansion 
bursts is more consistent with the phenomenol- 
ogy of the slow oscillators. However, conclusive 
assignment of the source to either class must wait 
until kHz QPOs are detected in the persistent 
emission. If SAX J1748. 9-2021 is a slow oscillator, 
then it has the highest frequency of any of the 
slow oscillators and the highest inferred spin fre- 
quency of any burst source, but it is slower than 
the spin frequency of the millisecond X-ray pul- 
sar XTE J1751-305 (Markwardt et al. 2002). If it 
is a fast oscillator, it has the lowest frequency of 
any of the fast oscillators and the lowest inferred 
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spin frequency of any burst source, but it is faster 
than the millisecond X-ray pulsar XTE J0929-314 
(Galloway et al. 2002). 

The dichotomy between slow and fast oscilla- 
tors may be related to the magnetic field configu- 
ration of the neutron star. Diametrically opposed 
magnetic poles provide a natural means to pro- 
duce oscillations with once or twice the spin fre- 
quency. However, the low upper limits placed on 
the subharmonic for fast burst oscillators are a se- 
rious problem for models containing two diamet- 
rically opposite hot spots because of the high de- 
gree of symmetry required (Muno et al. 2002) . The 
problem is particularly severe when oscillations are 
detected early in the burst rise because nearly si- 
multaneous ignition of the hot spots would be re- 
quired. 

The evolution of the magnetic field of a neutron 
star which is spinning-up or down is determined by 
interactions between superfluid neutrons and su- 
perconducting protons within the stellar core (An- 
derson & Itoh 1975). For a spinning-up star, quan- 
tized magnetic flux tubes in the superconducting 
core move inwards towards the rotation axis (Rud- 
crman 1991). For the NS-LMXBs of interest here, 
it is likely that sufficient accretion has occurred 
during spin-up so that the neutron star crust has 
been replaced several times, giving a correspond- 
ing movement of the surface magnetic fields even 
if neither eddy current dissipation nor creep in the 
crust has been effective. Since accretion is still oc- 
curring, the surface magnetic field emerging from 
the crust should closely follow the configuration of 
the core magnetic field. 

For the case where all of the emerging mag- 
netic flux returns to the same spin hemisphere, 
the magnetic flux is squeezed into a small polar 
cap where the spin axis intersects the crust (Chen 
et al. 1998). For spin-up from periods near 10 s to 
periods of milliseconds, both the north and south 
magnetic poles move close to one rotational pole 
and the net magnetic dipole field becomes nearly 
orthogonal to the rotation axis - an orthogonal 
rotator. The local surface magnetic field strength 
within the magnetic polar cap is then much higher 
than the dipole field, typically 10 11 - 10 12 G for 
a dipole field of 10 8 G. For the case where flux 
returns in the opposite hemisphere, the magnetic 
field is confined to small polar caps at the north 
and south spin poles. The dipole magnetic field in 



this spun-up configuration is oriented nearly par- 
allel to the spin axis - an aligned rotator. 

Radio millisecond pulsars (MSPs) are similar to 
NS-LMXBs, indeed they may be the descendants 
of the NS-LMXBs, in that they have experienced 
significant spin-up subsequent to their formation. 
The orientation of the magnetic field in radio pul- 
sars can be constrained using pulse profile and po- 
larization measurements. The fraction of MSPs in 
the Galactic disk which are either orthogonal ro- 
tators or nearly aligned rotators (exactly aligned 
rotators would not be pulsars, so there must be 
a slight offset, ~ 20°, of the magnetic field from 
the rotation axis) is significantly higher than for 
canonical pulsars (Xilouris et al. 1998). Chen et al. 
(1998) showed that this overabundance of orthog- 
onal and nearly aligned rotators in MSPs can be 
understood as a consequence of the magnetic field 
evolution of a neutron star with a superfluid and 
superconducting core, as described above. The or- 
thogonal and nearly aligned MSPs represent the 
two endpoints of the magnetic field evolution of a 
spun-up neutron star. 

The magnetic field configuration in LMXBs 
should be even more tightly compressed towards 
the rotation poles than in MSPs because in 
LMXBs the neutron star is still actively accret- 
ing and spinning up. The same physical processes 
which cause the migration of the magnetic poles 
toward the rotation axis during spin up cause 
migration away from the rotation axis during 
spin down. Because MSPs are spinning down, 
some motion of the magnetic field away from the 
poles has occurred (the magnitude of the motion 
is small, since the magnitude of the spin down is 
small compared to that of the spin up), thereby de- 
creasing the level of symmetry between the poles. 

It is interesting to speculate that the same di- 
chotomy may apply to NS-LMXBs. The orthogo- 
nal rotator geometry provides a natural means to 
produce oscillations at twice the spin frequency, 
while a slightly misaligned parallel rotator would 
produce oscillations at the spin frequency. The di- 
chotomy between the fast and slow burst oscilla- 
tors would then be related to physical differences 
between the magnetic field configurations of the 
neutron stars. 

The magnetic field geometry of orthogonal ro- 
tators would produce a high degree of symmetry 
between the two magnetic poles which would be 
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located very close to each other. This would natu- 
rally suppress the signal at the 1/2 subharmonic, 
i.e. the spin frequency, in the fast oscillators and 
allow nearly simultaneous ignition since the two 
poles are located very close together. Chen et al. 
(1993) estimate that spin up from 10 s to millisec- 
ond periods compresses the magnetic field into a 
region with a radius of about 10 4 cm around the 
spin axis. For a typical neutron star radius of 
10 6 cm, the poles would then be less than 1° from 
the spin axis and be antipodal within 2° as Muno 
et al. (2002) conclude is required to explain the 
lack of harmonics and sub-harmonics in burst os- 
cillations. The requirement that there be no more 
than 2% difference in the relative brightness of the 
two poles in the main and decaying portions of 
bursts (Muno et al. 2002) is more difficult to ad- 
dress, but the high degree of symmetry between 
the two poles and their close location for orthogo- 
nal rotators suggest that such uniformity may be 
achieved several rotation periods after ignition of 
the burst. 
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duty scientists of the BeppoSAX Science Opera- 
tions Center in the near to real-time WFC data 
analysis and the efforts of the RXTE team, par- 
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thanks Mai Ruderman for useful discussions and 
acknowledges partial support from NASA grant 
NAG5-7405. JZ acknowledges financial support 
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